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Abstract 

Front propagation described by Huygens' principle is a fundamental mechanism of spatial spreading of a property or an effect, 
occurring in optics, acoustics, ecology and combustion. If the local front speed varies randomly due to inhomogeneity or motion of 
the medium (as in turbulent premixed combustion), then the front wrinkles and its overall passage rate (turbulent burning velocity) 
increases. The calculation of this speedup is subtle because it involves the minimum-time propagation trajectory. Here we show 
mathematically that for a medium with weak isotropic random fluctuations, under mild conditions on its spatial structure, the 
speedup scales with the 4/3 power of the fluctuation amplitude. This result, which verifies a previous conjecture while clarifying its 
scope, is obtained by reducing the propagation problem to the inviscid Burgers equation with white-in-time forcing. Consequently, 
field-theoretic analyses of the Burgers equation have significant implications for fronts in random media, even beyond the weak- 
fiuctuation limit. 
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] 1. Introduction 

Phenomena from combustion [1] to seismic waves [2] to 

' population spreading [3] can be modeled by Huygens' prin- 

. ciple of front propagation, first stated as a law of geometri- 

' cal optics. In each application, the boundary of the affected 

. region (idealized as a sharp front) advances normal to itself 

■ at a locally specified speed. In a uniform medium, where 
, this speed is constant, an initially wrinkled front flattens 
' out over time; but in a spatially varying medium, a com- 
, petition occurs as wrinkling is continually reintroduced [4] . 

■ A central problem for the latter case is to determine the 
, overall statistically steady propagation rate of the wrinkled 
' front, which exceeds the average local speed because the 

front is defined by the fastest paths through the medium 
' (first passage). Our main result, establishing under general 
, conditions the proportionality of the speedup to the 4/3 
power of the amplitude of weak fluctuations, agrees with 
previous heuristic analysis [4] and numerical simulations [4- 
7], and eliminates ambiguities [8] concerning the physical 
relevance of the scaling. In the important case of premixed 
fluid combustion, the result describes the weak-advection 
scaling of the turbulent burning velocity [9] . 
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Our analysis relates the front dynamics to the evolu- 
tion of a pressure-free fluid obeying the white-noise-driven 
Burgers equation, itself a widely studied model of turbu- 
lence [10]. The results reported here enable adaptation of 
existing treatments of Burgers turbulence [11,12] to esti- 
mate the prefactor of the speedup scaling and its depen- 
dence on medium structure, with implications even for the 
opposite, practical limit of strongly advected flames. 

The propagation of a Huygens front in a general nonuni- 
form medium is governed by the local advecting velocity of 
the medium as a function of time and space, u(t, x) , and the 
local speed of propagation relative to the medium, v{t, x). 
Huygens' principle can then be stated as follows (using 
three-dimensional language for definiteness): If at time Iq 
a point xo lies in the affected region (including its bound- 
ary, the front), then at a slightly later time ti = + dt 
all points in a ball of radius w(to,xo)dt about the point 
xq + u{to, xq) dt are affected. The boundary of the affected 
region at ti (the new front) thus consists of certain points 
on the surfaces of balls originating from the initial front. 
If we consider all points on these spherical surfaces, then 
we automatically include those on the new front and more. 
Hence the front at any later time will be found among the 
affected points reached by arbitrary trajectories x(i) that 
start on the initial front and always move at the local speed 
V (in any direction) relative to the medium, so that they 
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obey 



- u(t,x) 



(1) 



Two simplifications of tliis general framework are phys- 
ically important. First, for geometrical optics in a static 
or "quenched" medium with nonuniform refractive index, 
or for combustion of a solid propellant with nonuniform 
burning rate, we have v = z;(x) (time-independent fluctu- 
ations) and u = (no advection). Note that our nonrel- 
ativistic equations do not correctly describe advection of 
light anyway, but are adequate for an optical medium "at 
rest" (where fluctuations of v dominate those of u) and 
for advection of sound in geometrical acoustics. Second, for 
idealized combustion of a premixed turbulent fluid in the 
limit of a very thin flame front [1], v is a constant (the lami- 
nar flame speed) and u(t, x) is the turbulent flow (assumed 
to be unaffected by the flame). Following the derivation of 
our key results in Sections 2 and 3, further implications for 
combustion are discussed in Section 4. 

2. Fronts and particles 

In the general case, Eq. (1) defines a large family of "vir- 
tual" trajectories x{t) of which only a subset actually form 
the front at a given time. The criterion for the relevant 
trajectories is simplest if |u| < v everywhere, as we now 
assume (in Section 4 we discuss relaxing this assumption). 
This inequality, which is trivially satisfied for a quenched 
medium, ensures that advection can never sweep the front 
backward and thus that each point y is crossed by the front 
only once. The time of this crossing, which we call To(y), 
is simply the time when the first virtual trajectory reaches 
y, hence defining a first-passage problem. 

To reduce the number of extraneous trajectories, we can 
exploit the first-passage criterion (minimization of travel 
time) and obtain constraints on relevant trajectories. Let 
us parametrize the solutions of Eq. (1) by 



dx 



u -|- vn, 



(2) 



with n(t) a unit vector. Then (sec Appendix A) a necessary 
condition for a first-passage trajectory is 



= -P„(A-n + Vt;), 



(3) 



where Pn denotes the projection orthogonal to n and A is 

the velocity gradient tensor Aij = Vi further neces- 
sary condition is that the trajectory starts out with n nor- 
mal to the initial front, from which it follows that n remains 
normal to the evolving front. In a quenched medium, where 
u and A are zero, Eqs. (2) and (3) reduce to the ray equa- 
tions of geometrical optics, and there it is well known that 
rays propagate normal to fronts. But with advection, we see 
that a trajectory's tangent vector dx/ dt is no longer aligned 
with the front normal n. As discussed in Section 4, Eqs. 
(2) and (3) govern front-tracking trajectories even when 




Fig. 1. Irrelevance of particles after collision. Two particles, starting 
from A and B on the initial front i, collide at C] the continuation of 
trajectory AC (dashed) reaches D. If it were the first arrival, its travel 
time would be an absolute minimum over all virtual trajectories from 
i to D, including BCD. Having the same minimum time, BCD must 
also obey the law of motion, but this deterministic law does not 
permit the time-reversed trajectory DC to split at C. Hence neither 
ACD nor BCD represents first passage to D, and the particles can 
be discarded upon colliding. 

|u| > f, but we assume |u| < v so that the front evolu- 
tion is described by a simple first-passage problem and a 
single- valued function To(y). 

For first-passage purposes, then, we consider a contin- 
uum of "particles" starting simultaneously from all points 
on the initial front (with initial n given by the unit normal) , 
and obeying Eqs. (2) and (3). If the Huygens front repre- 
sents a wave phenomenon in the geometrical-optics limit 
of very short wavelength, then we are describing physical 
quasiparticles: photons for light or phonons for sound. Mo- 
tivated by this physical case, we call Eqs. (2) and (3) the 
"law of motion" for first-passage trajectories. In premixed 
combustion, where the front represents a thin flame, the 
"particle" trajectories are mathematical constructs known 
as ignition curves [13]. 

Even with the condition (3), not all of these particles re- 
main on the front. The departure of particles from the front 
(into the interior of the affected region) is associated with 
"cusps" at which the normal n is not unique. Such cusps 
develop during propagation in a random medium even if 
the initial front is smooth [14,15]. When two distinct parti- 
cles reach the same point at the same time (necessarily with 
different n), both colliding particles fall behind the front 
and can then be discarded, as shown in Fig. 1. (In optics 
and acoustics, the corresponding photons and phonons re- 
main observable as second and later arrivals, but our scope 
is limited to first passage.) 

The collision rule implies a correspondence between our 
continuum of particles and a model "fluid" without pres- 
sure or viscosity, consisting of fluid elements that move in- 
dependently in response to external forces but disappear 
when they collide. Such a fluid is described by the inviscid 
Burgers equation (with suitable forcing) and the collisions 
are known as shocks [13,16]. Because our particles flll a sur- 
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Fig. 2. Complementary descriptions of first passage. A front propa- 
gates upward amid possibly time-dependent variations in u (arrows) 
and V (grayscale, dark for slow). We show u and v at the local time 
Tq. Snapshots of the front (thick curves) are contours of Tq; parti- 
cles (thin curves) track the front until they collide at cusps. Two 
circles expanding from affected points illustrate Huygens' principle. 
The fluctuations are weak enough to define a "fluid" using as 
"time" for the particles. For very weak fluctuations, the front re- 
mains nearly flat, and ux becomes irrelevant as it merely shifts the 
front without advancing it. 

face and not all of space, we must define the corresponding 
Burgers fluid more precisely. The particles, by definition, 
always represent the first arrival at their locations; thus if 
two of their paths reach the same point x, the particles 
necessarily arrive at the same time Tq (x) and are then dis- 
carded. Let us take the initial front to bo planar and use 
coordinates x = (a;||,xj_) such that this plane is a;|| = 0. 
Then the spatial paths of the particles can be viewed as 
"trajectories" x_L(a;||) in the "time" x\\. 

In this picture, a Burgers fluid exists in the lower- 
dimensional space , and the physical time information is 
carried by the fmiction ro(x|| , x_l). Here we must introduce 
the assumption that the medium fluctuations are weak, 
i.e., that |u| and the changes in v are both vanishingly 
small compared to the average value of v. In this limit, 
because trajectories cannot deviate significantly from the 
a;|| -direction without falling behind the front, the parti- 
cles necessarily collide before they are deflected enough to 
make the function xj_(x||) ill defined. 

The law of motion for these trajectories x_L(a;||) could 
be derived from Eqs. (2) and (3), but it is simpler to ob- 
tain the Burgers equation in its conventional Eulerian form. 
The same minimization principle that yields Eqs. (2) and 
(3) also implies that To(x) satisfies a generalized "eikonal" 
equation (see Appendix B) 

u(ro(x),x) • Vro(x) +t;(ro(x),x) |Vro(x)| = l. (4) 

The relation between Tq and the particle trajectories is il- 
lustrated in Fig. 2 for a two-dimensional medium. Because 
To is constant over a front, VTq lies along the front nor- 
mal n (in the direction of propagation), and so the physical 



velocity of particles in Eq. (2) can be written 
dx VTo 



(5) 



We next show that in the limit of weak random fluctuations, 
Eq. (4) reduces to the forced inviscid Burgers equation for 
a fluid consisting of these particles. 

3. Weak-fluctuation limit 

In a reference frame where the average value of u is zero, 
and in units such that the average value of v is unity, let us 
parametrize the weak fluctuations by 



u{t,x) = eU(ef,a;||,x_L), 
v{t, x) = 1 -|- e V{et, xu , xj_) 



(6) 
(7) 



where U and V are homogeneous isotropic random fields 

and e is taken asymptotically to zero. The time dependence 
of u and v is scaled by e because the natural source of 
time dependence is advection by u itself, which goes to zero 
with e. Heuristic scaling analysis [4,6] yields the following 
conclusions in the e ^ limit; The medium can be con- 
sidered effectively frozen (U and V time-independent); the 
advection component orthogonal to the overall prop- 
agation direction is irrelevant; the front reaches a statisti- 
cally steady state over a distance of order e~^/^; and the 
steady passage rate exceeds unity by an amount of order 

g4/3_ 

Guided by these expectations, we define rescaled quan- 
tities 



t($,x_l) = e~P[To{xu,x±) - Xu], 



(8) 
(9) 



where wc anticipate that p = | produces a useful e ^ 
limit, but we also consider slightly different p values to see 
why I is special. The point of the rescaling is that if V^r 
reaches a steady-state average value — C after a finite-^ 
transient, then V||To averages to 1 — Ce^^ after a charac- 
teristic distance xy e"^. A planar front propagating at 
constant speed would have VnTo = l/^-"*, and so to lead- 
ing order the passage rate in the random medium is 



(10) 



Our main result is a demonstration that for p = |, the 
rescaled front is governed by a forced Burgers equation that 
reaches such a steady state, thus implying e^/^ dependence 
of the speedup. 

Substituting Eqs. (6)-(9) into the eikonal equation (4), 
we find 

e C/|| (e^-f e + e^+f T, x^) (1 + e^^ V^r) 

+ [l + e Vie^-P^ + e^+Pr, e'P^, x^)] 

X ^'(l + e2pVjr)2 + e2p| V_lt|2 = 1. (11) 
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Fig. 3. Convergence of weak fluctuations to white noise. Choosing 
a simple random medium with (7(0)7(2)) ~ exp(— 2^), we show 
the two-point correlation function of rj, Eq. (13), for e = 1 (a) and 
€ = 0.2 (b). As € — > 0, the function is infinitely compressed in the 
"time" direction ^, but with the choice of scaling exponent p = ^, 
its time integral at a given z± (shaded area) remains fixed. This 
integral gives the zx-dependent coefBcient of a delta function that 
describes spatially correlated white-in-time noise. 

For i < p < 1, this result (multiplied by e"^^) simplifies 
considerably as e ^ 0: 

Vjr + i| V^rl^ = -ei-2f (0, e-^^ x^) 

= 77(e,xx). (12) 

As expected, both the component U_l and the time de- 
pendence of the medium go away in this limit. The front re- 
sponds to the frozen random field 7(.X||,xx) = ^{U\\ + V), 
which is homogeneous but not generally isotropic because 
a specific component of U is selected. The rescaled function 
??(?)X_l) = e'^"^'' 7(e~P^, xj_) appears to diverge as e ^ 0, 
but its dependence on the "time" ^ also becomes more and 
more rapid, and in fact with p = | it converges under mild 
assumptions to a white-noise process in ^. The two-point 



correlation function of rj is 

(r/(C, xx) r/(e + C, x± + zi)) = (7?(0, 0) r,(C, z^)) 
= e2-3Pe-f(7(0,O)7(e-^'C,zx)) 

/OC 
da(7(0,0)7(a,zx)), (13) 
-oo 

where the rescaled expression is squeezed into a Dirac delta 
function in ^, corresponding to white noise, times a factor 
that remains finite for p = |. Because any interval of ^ cor- 
responds to an infinitely long sample of the medium, r/ is a 
Gaussian random field even if 7 is not, provided the corre- 
lations of 7 decay fast enough with so that the central 
limit theorem is applicable. (For example, in homogeneous 
turbulence with strong intermittency, the statistics of 7 
are highly non-Gaussian; but because the flow dccorrclatcs 
over a few integral scales, the long-term effect on the front 
is given as e — > by Gaussian white noise.) Despite becom- 
ing white in "time" , the ry correlation function (13) remains 
smooth in "space"; this e ^ limit is visualized in Fig. 3. 

To confirm that Eq. (12) with p = | describes a Burgers 
fluid, note that 

VTo = (1 + e^/^V^T, £2/3 v^r), (14) 
and so the physical particle velocity (5) gives 
_ dx_L _2/3 dxj_/dt 
dx\\/dt 

-2/3 eU± + (l + ey)eV3v^T/|VTo| 
^ eC/|| + (l + ey)(l + e4/3V5r)/|VTo| 

^ VxT (15) 

for the Burgers "velocity" w in the time variable ^. The x_l 
gradient of Eq. (12) then yields 

V5W (w Vx)w = V_L7?(^, xx), (16) 

which is the inviscid Burgers equation with white-noise 
forcing [17] . Since the fluid elements are precisely the front- 
tracking particles, cusp singularities in the front are de- 
scribed by the standard jump conditions for Burgers shocks 
[18]. 

Because rj averages to zero, Eq. (12) implies that in a 
steady state 

C^-{V^r) = \{\V^T\^) = \{w^). (17) 

Thus the prefactor of e^/^ in the front speedup is the steady- 
state energy density of the white- noise-driven Burgers fluid. 
It was shown rigorously [17] that for white noise that is pe- 
riodic in space with any given smooth correlation function, 
there exists a unique statistically steady state for the invis- 
cid Burgers equation. This holds in any spatial dimension, 
generalizing a previous result [19] for the one-dimensional 
Burgers equation (corresponding to two-dimensional prop- 
agation) . Since an assumption of spatial periodicity is an ac- 
cepted device for treating bulk properties, we have convinc- 
ing evidence that C is well defined for a realistic smooth ran- 
dom medium given its two-point spatial correlation func- 
tion /(z) = (7(0) 7(z)), which determines the 77 correlation 
function (13), independent of higher statistics. 
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For geometrical optics in a quenched random medium, 
the essential rescaling technique and a reduction to white 
noise were previously demonstrated at the level of the ray 
equations [14,20], but the first-passage problem and the dis- 
appearance of rays at cusps (leading to the inviscid Burg- 
ers equation) were not considered. On the other hand, the 
analogous shock singularities of the Burgers fluid can be 
eliminated by a variant of Eq. (16) with an additional term 
— i^V^w on the left, known as the viscous Burgers equa- 
tion [11]. (The viscosity f > smooths the shocks.) The 
corresponding variant of Eq. (12), with — i/V^r on the left, 
is a form of the Kardar-Parisi-Zhang (KPZ) equation for 
interface growth [21]. Under forcing that is continuous in 
space and time, the solution of these equations is known 
to approach as ^ a limiting "viscosity solution" that 
reproduces the inviscid shocks [22] . Under white-noise forc- 
ing, the steady state of the viscous Burgers equation was 
analyzed using field theory [10,11], but the ly limit is 
more subtle. Fortunately, a recent theorem [23] establishes 
that the viscous white-noise steady state converges to the 
inviscid one, making the v ^ results of field theory ap- 
plicable to the e ^ limit of Huygens propagation. A pre- 
vious description of weakly advected flames based on the 
KPZ equation [9] assumed a white-noise velocity field for 
convenience, but did not observe that it could arise from a 
realistic smooth medium by rescaling, yielding the e*/^ law 
in the process. 

4. Conclusion 

The newly established equivalence between i/ — > Burg- 
ers turbulence and the e ^ propagation problem allows 
adaptation of field-theoretic [11] and numerical [12] meth- 
ods previously applied to the former in order to analyze 
the latter. Results of this adaptation will be reported 
elsewhere. The results obtained here already indicate the 
central importance of the spatial structure of the ran- 
dom medium to front propagation. Indeed, a sufficiently 
anisotropic medium can have a different scaling law, ex- 
plaining the dependence of the speedup derived for one 
such medium [8]. There, a flow U is constructed by adding 
(with random phases) a finite number of Fourier modes, 
none of whose wavevectors are exactly orthogonal to x\\. 
In this homogeneous but anisotropic medium, the integral 
in Eq. (13) oscillates and averages to zero. The amplitude 
of the white noise then vanishes and so C = 0, indicating 
correctly that the speedup goes to zero faster than e^/"*. 
As we will subsequently report, scaling in this random 
medium is a direct result of the gross anisotropy similar 
to the original argument for 6^ scaling [24], which was 
validated for an ccy -periodic medium [25]. We have here 
established e^/^ scaling as the generic case, applying to 
random media with a finite spectral density orthogonal 
to xy, such as isotropic media. (Although isotropic flows 
U lead to anisotropic 7, the generic scaling still holds be- 
cause in realistic cases the relevant transverse modes of 



have a finite spectral density; in particular, they are not 
constrained even if incompressibility is assumed. Further- 
more, generic quantitative deviations from flow isotropy 
that maintain this finiteness do not alter the scaling.) 

There is a different sense in which scaling exists even 
for isotropic media: as a transient starting from a flat ini- 
tial front. Before cusps form, a linearized approximation is 
adequate; we expect front tilts proportional to e and thus a 
speedup proportional to e^, both systematically increasing 
with time. In the Burgers picture, there is an order-unity 
interval of ^ before shocks become important, during which 
the steady energy input results in an energy density pro- 
portional to ^. The transient speedup is therefore propor- 
tional to e^/^^ = e^xii . As stated previously, our analysis 
applies after a distance x\\ ~ e~^/^, where a steady-state 
e*/^ speedup is attained. 

For premixed combustion applications, while weak ad- 
vection is a useful limiting case, the strong-advection limit 
|u| ^ t; is more important [1]. Because the white-noise re- 
duction and effective freezing of the medium no longer ap- 
ply, we expect the strongly turbulent burning velocity to de- 
pend not only on the two-point spatial correlation function 
/ but also on more complicated spatiotemporal flow statis- 
tics. A widely used isotropic-turbulence burning-velocity 
model [26] fails to incorporate this nonuniversality and also 
reduces to instead of e*/^ scaling in the weak limit. 

A further implication of our analysis for strong turbu- 
lence comes from a monotonicity property: For a given 
initial front and flow u(t,x), all fluid elements burned by 
a later time for laminar flame speed Vi are also burned 
for V2 > Vi. Thus the turbulent burning velocity can 
only increase with v, even into the weak-advection regime 
V ^ [u|. The aforementioned field-theoretic constraints on 
the weakly turbulent burning velocity then establish use- 
ful boundedness properties of the strongly turbulent one. 
Finally, although we derived the trajectory equations (2) 
and (3) assuming |u| < v, they are Galilean invariant and 
are thus valid in general, because Huygens' principle is lo- 
cal and we can switch to a reference frame where |u| < u in 
a neighborhood of a given event. The trajectory equations 
do not require a continuous representation of the front for 
constructing the normal n (except initially) and therefore 
may be useful in numerical studies at all turbulence inten- 
sities, although for very low intensity it is advantageous to 
exploit the Burgers equivalence. 
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Appendix A. Law of motion for first-pcissage 
trajectories 

The first arrival at an arbitrary point from an initial 
front is a constrained optimum: the result of minimizing 
the travel time T among virtual trajectories x(i) obeying 
Eq. (1), or equivalently (2) with arbitrary unit n(t), such 
that a = x(0) is on the initial front and y = x(T) is the 
desired endpoint. Using the method of Lagrange multipli- 
ers, we eliminate some of the constraints by introducing 
auxiliary variables. Consider a trajectory x(t) obeying the 
constraints. If the travel-time variation 5T vanishes for all 
infinitesimal trajectory variations (5x(f) that preserve the 
constraints (a necessary condition for a minimum), then for 
trajectory variations that keep a on the initial front but are 
otherwise arbitrary, 5T must be proportional to the devia- 
tions from the other constraints: 

8T = 5y-K+ I dt(5[u(i,x) + t;(i,x)n(i) -x(t)] • A(t). 
Jo 

(A.1) 

Here the overdot denotes a time derivative, k is the La- 
grange multiplier for the constraint on the endpoint y and 
\{t) is the Lagrange multiplier for the time-dependent con- 
straint (2). Equation (A.l) must hold for all infinitesimal 
variations (5x(t), (5n(f) as long as the starting point a 
remains on the initial front and n(t) remains a unit vector. 

The total variation of y = x(r) is Sy = 5^{T)+-k{T) 6T 
(i.e.. the endpoint can be changed cither by altering the 
trajectory itself or by following it for a longer or shorter 
time). Accordingly, Eq. (A.l) becomes 

ST = [(5x(T) + i(r) ST] ■ K 

+ I dt[Sx{t) ■ A{t,x) + Sx{t) ■Vv{t,x)n{t) 
Jo 

+ v{t,x)Sn{t)-S±{t)]-X{t). (A.2) 

Because ^n is constrained only by (5n • n = 0, the term 
proportional to (5n • A shows that A must be along n, say 
X{t) = n(t)/x(t). Next we integrate — i5x • A by parts to 

obtain 

ST = x(T) ■kST + Sx{T) ■ [k - n(T) /i(r)] 
+ SsL-n{0)fj,iO) 

+ [ dt Sx{t) ■ [A{t, x) • n{t) iJ.{t) + Vv{t, x) ij,{t) 
Jo 

+ h{t) n{t) + n{t) fi{t)]. (A.3) 

Validity of Eq. (A.3) for arbitrary ST gives 

1 = x(T) • K = [u(T, y) + v{T, y) n(T)] ■ k. (A.4) 

Validity for arbitrary Sx{T) gives 

K = n{T)fi{T). (A.5) 

Validity for arbitrary Sa. tangent to the initial front shows 
that n(0) is in the normal direction, as claimed. Because any 
instant of time can be considered "initial" , a first-passage 
trajectory always has n normal to the evolving front (while 



the trajectory remains on the front). Finally, validity for 
arbitrary Sx{t) gives (upon projecting orthogonal to n) 

n(i) = --Pn(t) (A(t, x) • n(i) + Vv{t, x)) , (A.6) 

which justifies Eq. (3) and, together with Eq. (2), defines 
the law of motion. Results equivalent to Eqs. (2) and (3) 
were obtained in geometrical acoustics by assuming, rather 
than proving, that n is the front normal [27]. 

Our variational method is reminiscent of the reciprocal 
Hamilton principle [28] , a modern contribution to the clas- 
sical mechanics of particles. Whereas the ordinary Hamil- 
ton least-action principle requires the action to be station- 
ary under variation of a trajectory between two points over 
a fixed time interval, the reciprocal Hamilton principle re- 
quires the travel time to be stationary under variation of 
a trajectory with fixed action. Both principles are equiva- 
lent to Newton's or Hamilton's equations of motion. In our 
problem, Eqs. (2) and (3) can in fact be derived from a 
Hamiltonian 

H{t,x,X) = u(t,x) ■ X + v{t,x.)\X\, (A.7) 

if we treat A (previously a Lagrange multiplier) as the 
canonical momentum and define n = A/jAj. The corre- 
sponding action is 

S = [ dt[X{t)-x{t)-H{t,x,X)] 
Jo 

= dt[±{t)-u{t,yi)-vit,x)n{t)]-X{t), (A.8) 

^0 

which vanishes for "Huygens" trajectories that obey Eq. 
(2), as well as for many other trajectories that do not. 

The reciprocal Hamilton principle implies that a first- 
passage trajectory has stationary travel time not only 
among Huygens trajectories but also among the wider class 
with 5 = 0. One might suppose we could consider the tra- 
jectory with absolute minimum (hence stationary) travel 
time in the 5 = class, apply the reciprocal Hamilton 
principle to obtain the equations of motion (2) and (3) and 
then observe that such a trajectory obeys the fundamental 
constraint (1) of Huygens propagation and thus represents 
the physical first passage. But this shortcut fails because of 
an important technicality: Among trajectories with 5 = 0, 
there is no absolute minimum travel time between two 
points. For simplicity, take a quenched medium with u = 0. 
A sufficient condition for 5 = is that x • A = v\X\, i.e., 
the component of x along A equals v. A suitable A exists as 
long as |x| > V, and so the 5 = class includes arbitrarily 
fast trajectories. Therefore the reciprocal Hamilton prin- 
ciple is not a sound basis for deriving the first-passage law 
of motion, although that law can be expressed in Hamil- 
tonian form. Instead, we have derived the law of motion 
from the correctly constrained variational principle. 

Appendix B. Generalized eikonal equation 

The extra implications of our variational principle give a 
useful result for the spatial dependence of the first-passage 
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time To{y). Equation (A.l) shows that for variations re- 
stricted to first-passage trajectories (T = Tq), so that the 
variation inside the integral is zero, we have STq = k ■ Sy 
and thus VTq = k. Equation (A. 5) and |n| = 1 then give 

VTo 



n 



I VTo I 



(B.l) 



confirming that n is always normal to the front. Substitut- 
ing into Eq. (A. 4), we find 

u(To,y) ■ VTo + «(To,y) |VTo| = 1. (B.2) 

A form of this equation is known in geometrical acoustics 
[29] as a generalization of the "eikonal" equation of geomet- 
rical optics. We see that it follows directly from Huygens' 
principle in a time-dependent advected medium. 

In fact, a closely related "G equation" was obtained for 
idealized premixed combustion [15,30]: 

G + u-VG = v\VG\, (B.3) 

where a particular "level surface" of the function G, say the 
set of points y where G{t, y) = 0, represents the front at 
time t. The gradient of the identity G(To(y),y) = gives 
VTo = —VG/G and proves the equivalence with Eq. (B.2). 
The G equation corresponds to the Hamilton-Jacobi equa- 
tion of classical mechanics [31], with Hamiltonian (A. 7), 
momentum A = — VG and Hamilton's principal function 
— G. Although the equations of motion for this Hamiltonian 
reproduce (2) and (3), in general such trajectories ("char- 
acteristics" of the Hamilton-Jacobi equation) do not move 
with level surfaces of Hamilton's principal function. They 
do so for this Hamiltonian because the action (A. 8) van- 
ishes, reflecting the equality of group and phase velocity 
(absence of dispersion) in Huygens propagation. 
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